Wave-unlocking transition in resonantly coupled complex Ginzburg-Landau equations 
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We study the effect of spatial frequency-forcing on standing-wave solutions of coupled complex 
Ginzburg-Landau equations. The model considered describes several situations of nonlinear counter- 
propagating waves and also of the dynamics of polarized light waves. We show that forcing introduces 
spatial modulations on standing waves which remain frequency locked with a forcing-independent 
frequency. For forcing above a threshold the modulated standing waves unlock, bifurcating into a 
temporally periodic state. Below the threshold the system presents a kind of excitability. 
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Different physico-chemical systems driven out of equi- 
librium may undergo Hopf bifurcations leading to rich 
spatio-temporal behavior. When these bifurcations oc- 
cur with broken spatial symmetries, they induce the for- 
mation of wave patterns described by order parameters 
of the form : 
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where the slow dynamics of the wave amplitudes A and 
B obey complex Ginzburg-Landau equations. This is the 
case, for example, for Rayleigh-Benard convection in bi- 
nary fluids, Taylor-Couette instabilities between corotat- 
ing cylinders, electro-convection in nematic liquid crys- 
tals Q , or for the transverse field of high Fresnel num- 
ber lasers . Symmetry breaking transitions are usually 
very sensitive to small perturbations or external fields. 
For example, it has been shown that a spatial modu- 
lation of the static electro-hydrodynamic instability of 
nematic liquid crystals modifies the selection and stabil- 
ity of the resulting roll patterns. In particular, the con- 
straint imposed by a periodic modulation of the instabil- 
ity point may lead to a commensurate-incommensurate 
phase transition In the case of Hopf bifurcations, 
external fields inducing spatial or temporal modulations 
strongly affect the selection and stability of the result- 
ing spatio-temporal patterns. For example, standing 
waves may be stabilized by purely temporal modulations 
at twice the critical frequency |^,^, or by purely spa- 
tial modulations at twice the critical wavenumber in 
regimes where they are otherwise unstable, including do- 
mains where the bifurcation parameter is below the crit- 
ical one. 

External forcings that break space or time transla- 
tional invariance, but not the space inversion symme- 
try of the waves amplitudes, induce linear resonant cou- 
plings between the complex Ginzburg-Landau equations 
(CGLE) which describe the dynamics of the amplitudes 
of left and right travelling waves. In the case of forcings 
that break the space translation invariance, the coupling 
coefficient e is in general complex, and the correspond- 
ing coupled CGLE may be written, in one-dimensional 



geometries, as : 

A + VgdxA ^ iiA+{l + ia)dlA 

B - Vgd^B = /iS + (1 + ia)dlB 

-{l + ^P){\B\^+^\A\^)B + eA (2) 

Due to the resonant coupling with coefficient e, pure 
travelling waves are not solutions any more of these equa- 
tions, and generic arguments of bifurcation theory allows 
a characterization of the possible uniform amplitude so- 
lutions depending on the various dynamical parameters 
of the system Here also, standing waves may be sta- 
bilized as the result of phase locking between the waves A 
and B. Predictions based on (^) in the a;— independent 
case have been successfully tested for azimuthal waves 
in an annulus laser with imperfect 0(2) symmetry [Q. 
However, the combined effect of the complex coupling 
coefficient e and the spatial degrees of freedom has not 
been explored. 

In this Letter, we study equations (2) with the fol- 
lowing parameter restrictions: imaginary linear coupling 
coefficient (e = i^p), negligible group velocity Vg, and 
weak and real nonlinear cross-coupling term (7 < 1). 
We will however maintain the spatial derivative in the 
r.h.s. of (HI), and this will be crucial for the results below. 
We will show that the spatial forcing introduces spatial 
modulations of the standing waves solutions while A and 
B remain frequency locked with a forcing-independent 
frequency. By increasing the forcing, these stable mod- 
ulated waves merge with unstable ones in saddle-node 
bifurcations with nontrivial global structure. This wave- 
unlocking transition results in a mixed state with limit 
cycle temporal behavior. The threshold value of the forc- 
ing and the limit cycle frequency are calculated analyti- 
cally. Modulated standing waves can also be induced by 
strong enough temporal forcing |^ . 

The parameter regime explored here would be appro- 
priate in physical situations where a spatial forcing mod- 
ulates the frequency of the Hopf instability and induces a 
purely imaginary resonant forcing (a purely real e would 
appear due to a spatial modulation of the distance to the 
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instability point). Possible systems should have negligi- 
ble group velocities, as in some circumstances in binary 
fluid convection ^ or liquid crystals ; and weak cou- 
pling such as in viscoelastic convection Up to now, 
the parameter range considered here best applies to sev- 
eral situations in laser physics. A first one corresponds 
to taking into account transverse effects in inhomoge- 
neously broadened (7 < 1) bidirectional ring lasers 
The purely imaginary resonant coupling is a consequence 
of conservative (off-phase) backscattering or alter- 
natively, a spatial modulation of the refraction index of 
the laser medium. In fact, a spatially periodic refractive 
index is the mechanism used for single frequency selection 
in index coupled Distributed Feedback lasers (DFB). A 
second situation is that of the transverse vector field in a 
laser near threshold The parameter 7^ corresponds 
to a detuning splitting between light linearly polarized 
in different orthogonal directions, produced for example 
by small cavity anisotropies. In this case, A and B are 
not the amplitudes of left or right travelling waves, but 
the amplitudes of the two independent circulary polar- 
ized components of light, that is A = [A^ + iAy)/ \pl and 
B = {Ax — iAy) 1^/2^ where A^ and Ay are the linearly 
polarized complex amplitudes of the vector electric field 
with a spatially transverse dependence. Weak coupling 
(7 < 1) favors linear polarization {\A\ = \B\). We will 
often use the light-polarization terminology, because it 
gives a clear physical insight into the states found for the 
general set of Eqs. (Q) of broad applicability within the 
parameter restrictions above. 

Two families of solutions of the coupled CGLE (|^) can 
be distinguished. The first family corresponds to trav- 
elling waves for A and B with the same amplitude, fre- 
quency, and wavenumber: 



(3) 



Without forcing (7P = 0), the constant global and rel- 
ative phases, 6*0 and ■0o, are arbitrary, the amplitude is 
Qg = (/i — fc^)/(l + 7), and the frequency w is = 
—ak"^ ~ + 7)Q§. With forcing, the global phase and 
the amplitude remain unchanged, but the relative phase 
is fixed by sin2'(/'o = 0; the two allowed values of V'o give 
two solutions with frequencies a; = wg ± 7p . The phase 
instabilities of these solutions were discussed in |l4]. 

The second family of solutions can be searched in the 
form of two waves 



A 



B = 



(4) 



frequency-locked to a frequency loq independent of forc- 
ing. For 7p = 0, the exact solutions of (|^) in this form 
only have two terms: |ai| — |5_i| = Qg- The effect 
of a small forcing in this solution is to generate higher 
harmonics, while keeping cjg fixed and the relative phase 



between ai and 6_i arbitrary. Now, the remaining co- 
efficients a„ and 6„ are not zero and can be calculated 
perturbatively in 7^. Close enough to the threshold for 
a mode k — k"^ « 0), the amplitude of higher order 
harmonics is negligible and, to lowest order in /i — /c^, an 
approximate solution takes the form 

B = e^i^o+^o^) [Qe~^^^^'^^°^ + Re''^''^'^^'''^'^'^ ) (5) 

with 60 and ^g arbitrary, and fixed by the forcing. 
Q and R are real numbers (positive or negative) and, 
for small 7p, \R\ <C \Q\ (an equivalent solution is found 
interchanging Q and R). 

A visualization of these solutions can be given within 
the polarization interpretation of (||) . Defining C±e*^± = 
Q ± e"^i?, the change of variables to the amplitudes of 
the X— and ?/— linearly polarized components gives 



Ax = ^/2C+ cos(fca: + ^g)e'('^ot+eo+f+) ^ 
Ay = V2C- sm{kx + V'g)e'('^»*+''«+f-) , 



(6) 
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FIG. 1. Top: Modulus (left) and phase (right) of A for 
7p — 0.012 < 7pc. The horizontal axis is space (256 units), 
and the vertical is time (1000 units). Gray levels range from 
black (0) to white (the maximum of the modulus or 2n for 
the phase). This numerical solution has been obtained from 
(2) with (3 = 0.2, 7 = 0.5, fj, = 0.2, a = 2.6. The initial 
condition is a standing wave with ui — uo and k = 0.123. 
Bottom: Polarization representation of the solution at a given 
point X. For 7p = one has linear polarization (indicated by 
the straight line) which becomes elliptical for 7p 7^ 0. An 
equivalent solution has the major axis of the ellipse along the 
second and fourth quadrants. 
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FIG. 2. Same as Fig. 1 but using 7p = 0.0145 > -ypc- 

These equations describe at each point x the superpo- 
sition of two dephased harmonic motions with different 
amplitudes and a frequency loq independent of forcing. 
This identifies the solution (H) with an elliptically polar- 
ized standing wave pattern in which the orientation of the 
ellipse and its cUipticity vary periodically in the spatial 
coordinate x. In the limit of no forcing, R = 0, the ellipse 
degenerates in a linearly polarized standing wave with an 
angle of polarization -0 — kx + ^Q. In this interpretation, 
the first family of solutions, (||), would correspond to lin- 
early polarized traveling waves with frequency w and an 
angle of polarization ?/;o- In such a case, the forcing fixes 
the direction of polarization so that only x- or j/-linearly 
polarized waves remain. On the contrary, forcing in (|^) 
grows an ellipse from a linearly polarized standing wave 
keeping the frequency unchanged. 

Elliptically polarized standing wave patterns are ob- 
tained from a direct numerical integration of the cou- 
pled CGLE as shown in Fig. 1. Increasing the forcing, 
these solutions become unstable through a bifurcation 
in which Q and R become time dependent. As shown 
in Fig. 2, the solution beyond this instability oscillates 
between the two equivalent elliptically polarized stand- 
ing wave patterns found for small 7p. In addition, from 
the numerical simulations, one finds that the period T of 
these oscillations decreases beyond the critical value 7pc. 
One has T"^ oc (7p — 7pc) (see Figs. 3 and 4). 

A quantitative description of the instability, includ- 
ing the determination of the critical forcing 7Pc and the 
period of the oscillations, can be performed by an am- 
plitude analysis. Close to the threshold for the k modes, 
the equations for the slow time evolution of Q and R 
can be found by substitution of into (||) and neglect- 
ing contributions from higher order harmonics. Defining 



X = {ii- k^)X - (1 +7)^3 - (1 + 7COS20)— sin2 2$ 



X^ 

$ = - (1 + 7 COS 2(?!)) — sin 2$ cos 2$ 
X'^ 

+ l3j sin 2'/'-^ sin 2$ — 7p sin ip , 

(/) = /3(1 + 7 cos 2(j))X'^ cos 2$ + 7^^ sin 2(/) 
— 2jp cos (f) cot 2$ . 



(7) 



The fixed points of (|7|) represent the polarized stand- 
ing waves solutions (P)T These points can be determined 
exactly in the limiting case of P — 0. The interesting so- 
lutions have two allowed values of (po = (2n + l)7r/2, 
n = 0, 1; and for each value of (j), there are eight fixed 
points: Four are stable {+) and the other four are saddle 
points (— ). The corresponding values of X and $ are: 



2(3 + 7)(l+7) 
5 + 3j± 



'(I-7F 



(l+7)(3 + 7)7j> 



(8) 



<i>o(+) = + M-) = 5 - e(-) + m| , (9) 



where m — 0,1, 2, 3, and 



(-1) 



ji+i 



47P 



(l-7)Xo(±)2 



(10) 



Heteroclinic orbits connect the saddles and the stable 
nodes with the same (p. When 7p grows, saddles and 
nodes approach by pairs and at the critical value 



7Pc = 



(M-fc^)(l-7) 
v/8(l + 7)(3 + 7) ' 



(11) 



they merge and disappear via inverse saddle node bifur- 
cations. The interesting point is the global structure of 
the bifurcation: the presence of the heteroclinic connec- 
tions gives rise to the birth of limit cycles (one for each 
value of 0). This is similar to the Andronov-van-der-Pol 
bifurcation |]l5|| that appears in several types of excitable 
systems [|l^. The difference is that, due to symmetries, 
here we have several pairs of fixed points merging, instead 
of just one pair. The periodic behaviour is illustrated 
in Fig. 2 by the periodic alternation of the trajectory 
between the "ghosts" of the disappeared elliptically po- 
larized states corresponding to the fixed points. Below 
the bifurcation, small perturbations around the stable 
solutions decay, whereas perturbations above a threshold 
push the system along the heteroclinic trajectory towards 
another stable fixed point. Since the size of the pertur- 
bation required for such switches decreases by increasing 
7p, and vanishes at "fpc, the multistability of this sys- 
tem can be seen as a kind of excitability |T7|] . A different 
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consequence of the multiplicity of stable states is their 
possible coexistence in space, leading to the formation of 
domains with different polarizations along the x axis. 




r^-cr^ 0.0140 0,0145, O.OIqQ 0,0155, , , „ , 

FIG. 3. Period or the osculations (T) of Q and it obtained 
from the numerical sdlution of the coupled CGLE using the 
parameters given in Fig. 1. The dashed line is a least squares 
fitting from which 7pc = 0.0138. 
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FIG. 4. The amplitudes of Q and R of the solution shown 
in Fig. 2 exhibit the periodic oscillations given by Eq. (13). 
The time has been scaled using the value of 7Pc obtained 
from Fig. 3. The dotted line corresponds to 7p = 0.0155, the 
dashed line to 7p = 0.0140, and the lines in between to the 
other points of Fig. 3. 

Close to the instability at 7Pc, the time dependent be- 
havior of the solution can be obtained reducing the prob- 
lem to a phase dynamics by elimination of the variable 
X. We have (in the limit /3 = 0), 




(/i- fc2)(l-7) sin 4$ 



IP 



(12) 



5 + 37 - (1 - 7) cos 4$ 
which for 7p > 7Pc yields the following time behavior: 

tan(2$) = tan(2$,) (l + v/2(7p - 7Pc)/7Pc 
(5 + 37)^7^ 



previous analysis is still meaningful as explicitly seen in 
the numerical results of Figs. 3 and 4. 

In summary, in the absence of forcing, and for the pa- 
rameter regime considered here, there are solutions for 
the amplitudes A and of the coupled CGLE which cor- 
respond to linearly polarized standing waves. We have 
shown that an imaginary coupling between them trans- 
forms these solutions into standing waves with spatially 
periodic elliptic polarization. Increasing the forcing, an 
instability of these solutions, via the unlocking of the un- 
derlying wave amplitudes, appears and the solutions ac- 
quire a time-periodic behavior. Locally, this bifurcation 
is of the saddle-node type, but the presence of hetero- 
clinic connections between the fixed points gives rise to 
the appearance of a limit cycle when stable and unstable 
points merge. 
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tan 



7P - 7Pc 



(l+7)(3 + 7) 



(13) 



where, $c = Eq. (|), for 7p = 7pc. 

An approximative analysis of Eqs. (|^) for /3 7^ in- 
dicates that this parameter appears squared in the ex- 
pressions for 7pc, X and <i>c. Therefore, for small /?, the 
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